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Using the partially completed Submillimeter Array with five antennas, we 
have observed the CO J=2-l and 3-2 emission from the envelope surrounding 
the carbon star V Hya. The high-angular resolution (2"-4") maps show that V 
Hya is powering a bipolar molecular jet having an extreme velocity of 70-185 
km s~^. The axis of this high velocity jet is perpendicular to the major axis of 
! the fiattened disk-like envelope, which is expanding with a velocity of ~16 km 

'. s~^. There is a third kinematic component, a medium- velocity wind having a 

deprojected velocity of 40-120 km s~^ moving along the disk plane. Both the 



O ! high velocity jet and the medium velocity wind have a dynamical time scale of 

! a few hundred years. The fiattened structure and the coUimated jet observed in 

r-| ! V Hya suggests that the formation of asymmetrical structure proceeds while the 

^! central star is still in the AGB phase. 

O 

H , 

! Subject headings: ISM: individual (V Hya) — ISM: jets and outfiows — ISM: 
molecules — stars: AGB and post-AGB 



."^ , 1. INTRODUCTION 

V Hya, an evolved star located at 380 pc away from the sun (Knapp et al. 1997), is 
a carbon star with a high mass loss. Although the optical properties of V Hya are those of 
the normal N type carbon star, the envelope surrounding this star shows several peculiar 
properties. Single-dish molecular line observations have shown that the envelope of V Hya 
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is elongated along the north-south direction with a velocity gradient along its minor axis 
(e.g., Kahane et al. 1996). Such morphological and kinematical properties are different 
from those of the typical asymptotic giant branch (AGB) envelope, which has spherically 
symmetric shape and uniform expansion. In addition, V Hya is associated with a very fast 
(> 100 km s^^) outflow that was observed in the optical and infrared spectra (e.g., Sahai 
& Wannier 1988; Lloyd Evans 1991), and in the CO J=2-l and 3-2 spectra (Knapp et 
al. 1997). Knapp et al. (1997) also found that the high- velocity CO gas is bipolar, and 
proposed that the fast-moving gas is expanding along the east-west direction, perpendicular 
to the major axis of the envelope. However, previous observations did not bring us the 
detailed structure of the high-velocity outflow and the relation between the high-velocity 
flow and the slowly expanding envelope. 



2. OBSERVATIONS 

The observations were carried out between February 2003 and May 2003, with the five 
antennas of the partially completed SubmiUimeter Array (SMA) at Mauna Kea, Hawaii 
(Ho, Moran & Lo 2004). The primary-beam size (HPBW) of the 6 m antennas at 230 and 
345 GHz were measured to be ~54" and ~36", respectively We obtained the data in two 
array configurations, which provided 15 independent baselines. The spectral correlator had 
a bandwidth of 656 MHz and a frequency resolution of 812.5 kHz. 

The visibility data were calibrated using the MIR package (Scoville et al. 1993). We 
used 1055-1-018 or 3C273 as a phase and amplitude calibrator, and Callisto as a flux cali- 
brator. The bandpass was calibrated by observations of Jupiter or Saturn. The calibrated 
visibility data were Fourier transformed and CLEANed by using the Astronomical Image 
Processing System (AIPS) with natural weighting. The synthesized beam had a size of 
4.4"x3.1" with a position angle of -21° at 230 GHz and 2.1"xl.9" with a position angle of 
—28° at 345 GHz. In making maps, we averaged 2 channels in CO J=2-l and 3 channels 
in CO J=3-2 and made maps with velocity resolution of 2.1 km s~^. The rms noise level of 
the maps at the velocity resolution of 2.1 km s~^ was 0.09 Jy beam~^ (0.15 K in Tb) in CO 
J=2-l and 0.3 Jy bcam^^ (0.77 K in Tb) in CO J=3-2. We averaged the entire 656 MHz 
of both image sidebands to make the maps of 241 GHz and 335 GHz continuum emission. 
Continuum maps were made with uniform weighting, which provided a synthesized beam of 
3.0"x2.2" with a position angle of 37° at 241 GHz and 2.0"xl.3" with a position angle of 
48° at 335 GHz. The rms noise levels of the 241 GHz and 335 GHz continuum maps were 9 
mJy beam~^ and 20 mJy beam~^, respectively. 
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3. RESULTS AND DISCUSSION 

3.1. Continuum emission from V Hya 

We detected the continuum emission peaks at q;j2ooo — 10*^ 51™ 37.25^, Sj2ooo — — 21°15' 
00.5" at both 241 and 335 GHz. The position of the continuum source coincides well with the 
stellar position at the observational epoch (J2003.25), which is derived from the coordinates 
and proper motion measured by the Hipparcos satellite (Ferryman et al. 1997). Figures Ik 
and 11 show that the continuum source is not resolved with our beams. This suggests that the 
most of the emission comes from a spatially compact region with a radius of <1" (380 AU). 
The flux density of the source is 63±16 mJy at 241 GHz and 105±39 mJy at 335 GHz. The 
spectral index a between 241 GHz and 335 GHz is ~1.6, which is shallower than Rayleigh- 
Jeans. This might imply the contribution of the free-free emission from shock-ionized gas 
close to the star (Luttermoser & Brown 1992; Knapp et al. 1997). However, because of the 
large uncertainty in the flux at 335 GHz, it is diflicult to estimate the contribution of the 
ionized component to the submillimeter spectra. If we assume that the continuum emission 
at 241 and 335 GHz arises from the thermal emission from the dusty envelope, the mass of 
dust Mdust required to produce the observed flux densities is estimated by 

n/T ^'2'apd,ustDX'^{S„ — S*) , , 

Mdust^ on i,^ ' 

in the Ray leigh- Jeans limit. Here, Si, is the observed flux, 5"* is the flux that comes from the 
stellar photosphere, and D is the source distance. At 241 GHz, is estimated to be ~10 
mJy We assume that the dust particles have radius a of 2000 A, material density Pdust of 
2.25 g cm^'^, and cmissivity of 6.17x10^"^ at 274.6 GHz with a power-law dependence on 
frequency with index [3=1. We adopted a grain temperature to be 360 K. This corresponds 
to the grain temperature at r =380 AU if it has a power-law dependence on radius, Tg = 
T!,(-R*/r)^/*^^+'^) (Knapp, Sandell, & Robson 1993), with an effective temperature and a 
radius i?* of the star of 2650 K and 3.8x10^^ cm (Knapp et al. 1999), respectively. The 
mass of dust was calculated to be 6.5x10"^ Mq. If we assume a constant mass-loss, the dust 
mass loss rate is estimated to be 5.7xlO~^M0 yr~^ based on an outflow speed of ~16 km 
s"-*^ (see section 3.4). 



3.2. Flattened envelope and high-velocity bipolar outflow 

Strong CO J=2-l and J=3-2 emission was detected in the velocity range from Vlsr = 
—45 km s~^ to -1-12 km s~^, corresponding to ±~28 km s~^ from the systemic velocity, 
Vgys =—17.5 km s~^. Both CO J=2-l and 3-2 maps show that the molecular envelope is 
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elongated along the north-south direction. We have smoothed our interferometric maps to 
mimic the CO J— 2-1 and 3-2 spectra observed with the 30" beam and 20" beam, respec- 
tively, and compared the spectra observed by Knapp ct al. (1997) and Stanek et al. (1995). 
We found that ~50 % of the CO J=2-l flux and ~35 % of the CO J=3-2 flux observed by 
the single-dish telescope was recovered by the SMA. 

The velocity-channel maps of CO J=3-2 emission are shown in Figures la to li. The 
velocity structure traced by the CO J=2-l emission is similar to that shown in the CO 
J=3-2 maps. We found that the approaching part of the envelope is located to the west and 
the receding part to the cast of the star, which is consistent with the previous single-dish CO 
J=2-l results of Kahane et al. (1996). The spatial- kinematic structure shown in Figure 
1 can be explained if the CO emission arises from an inclined flattened envelope expanding 
along its radial direction. A flattened disk-hke structure was also identified by Sahai et 
al. (2003) based on their interferometric CO J=l-0 observations. However, the observed 
velocity pattern does not completely match that of a uniformly expanding disk. The largest 
spatial extents of ±7" toward the east and west are shown in the velocity channels of ±~8 
km s~^ from Vgys (Figures Ig and Ic, respectively), and not in the highest velocity channels. 
On the other hand, the higher velocity CO emission shown in Figures la and li comes from 
the regions close to the central star. The position-velocity diagram of CO J=2-l along the 
east-west cut presented in Figure 2b suggests that there are two kinematic components, that 
is, two vertical ridges at Ay~±8 km s~^ from Vgys and a horizontal ridge that extends to 
±30 km s"^. The ±8 km s""^ components correspond to the two spikes in the CO hue profiles 
observed by Kahane et al. (1996) and Knapp et al. (1997), and are hkely to arise from the 
expanding disk inclined to the line of sight (hereafter referred to as the low- velocity disk). 
The ±30 km s~^ component has a clear bipolarity with the redshiftcd part to the east and 
the blueshifted part to the west (referred to as the medium- velocity wind). 

Both CO J=2-l and 3-2 lines are associated with the broad wing emission that extends 
to ±~150 km s~^ from Vgys- The detection of these wing components is confident because 
the phases of the visibilities show coherent behavior across the velocity range of ±150 km 
s~^. In Figures 2a and 2c, we show the position- velocity diagrams of the CO J=2-l along 
the east-west cut in the blueshifted {AV from Vgyg < —40 km s~^) and redshifted (AV^ > 
±40 km s~^) velocity ranges, respectively. In order to improve the signal to noise ratio, we 
made the maps of these velocity ranges by averaging 10 channels, providing the effective 
velocity resolution of 10.6 km s~^. In order to exhibit the features clearly, we smoothed the 
maps to 5" before we create these three position-velocity diagrams. 

Figures 2a and 2c suggest that emission in the velocity ranges of AV < ±60 km s~^ 
corresponds to the lower intensity part of the medium-velocity wind. On the other hand. 
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the velocity structure of the higher velocity component with AV > ±60 km is different. 
The spatial distribution of the high-velocity component {AV from V^ys is ±(60-162) km 
s~^) shown in Figure 3 superposed on the intensity- weighted mean velocity map of CO 
J=2-l. The high- velocity component has its approaching part in the east and the receding 
component in the west of the star. This orientation is opposite with respect to the lower 
velocity components as suggested by Knapp et al. (1997). Such an orientation can be 
explained if the high velocity component is a pair of coUimated jets blown out in the polar 
direction of the low-velocity disk. Hereafter, we refer to this component as the high- velocity 
jet. Recently, Sahai et al. (2003) discovered the [S II] emission-line blob with high negative 
radial velocity of Ay~— 240 km s~^ at 0.2" east of the central star. This also suggests the 
presence of the collimatcd high velocity jet with approaching part to the east of the central 
star. In the CO J=2-l, the receding part of the bipolar jet, which has not observed in 
the [S 11] emission probably due to the absorption caused by the obscuring disk, is clearly 
detected in the west of the star. Figures 2a and 2c show that the negative and positive 
velocities of the high velocity jet increase as the positional offset from the star increases. 
Such a velocity increase can be explained if we take into account that the density of the 
circumstellar envelope decreases rapidly along the east-west direction. Figures 2a, 2c and 
3b show that the high velocity component consists of discrete clumps separated by ~5", 
suggesting that the high velocity jet is episodic. 



3.3. Circumstellar structure of V Hya 

The high-angular resolution CO J=3-2 and 2-1 maps suggest that the circumstellar 
structure of V Hya consists of three kinematic components; the low-velocity disk having 
a radial velocity offsets of AV — ±8 km s~^ from Vgys , the medium-velocity wind with 

AV < ±60 km s~^, and the high- velocity jet with AV^=±60-160 km s~^. These three 
components corresponds to the ±8 km s~^ spikes, the 45 km s~^ wind, and the 200 km s~^ 
wind defined by Knapp et al. (1997), respectively. 

The spatial extent of the low-velocity disk along the major (N-S) and minor (E-W) axes 
measured at the 3 a level are ±8" and ±7", respectively in the CO J=3-2 and ±15" and 
±13", respectively in the CO J=2-l. If we assume that the disk is thin, the ratio between the 
major and minor axes suggests that the disk plane is inclined by ~60° from the line of sight, 
which is consistent with the the previous estimation of ~70° by Sahai et al. (2003). We find 
that the gas in the low-velocity disk is moving at almost constant deprojected velocity of 16 
km s~^, which is similar to the typical AGB wind velocity. The north-south extent in the 
CO J=2-l map suggests that the low- velocity disk has a radius of 5700 AU and a dynamical 
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time scale of 1700 yr. If we assume that the axis of the high velocity jet is perpendicular to 
the disk plane, the deprojected velocity and the dynamical time scale of the jet are estimated 
to be 70-185 km s~^ and 100-250 yr, respectively. 

What is the nature of the medium- velocity wind? Here we present three possible scenar- 
ios; 1) the interaction between the low-velocity disk and the high-velocity wind, 2) equatorial 
mass outflow with higher velocity, or 3) a second bipolar outflow with different axis. The first 
case assumes that the dense equatorial disk collimates the spherically expanding or poorly 
collimated high-velocity wind into the polar direction. In this case, the wind directed to 
the polar direction can escape rather freely, while the one directed to the equatorial plane 
pushes the dense molecular gas in the disk plane, and produces the medium-velocity wind. 
The second case is possible if the velocity of the recent mass fiow along the equatorial disk 
plane has been remarkably enhanced for some reason. In either case 1 or case 2, the di- 
rection of the medium-velocity wind is close to the disk plane. Therefore, the deprojected 
velocity of the wind is estimated to be 40-120 km s~^. The dynamical time scale of the 
medium-velocity wind estimated from the spatial extent of ~5-8", inclination angle, and 
the deprojected velocity is 100-400 yr, which is close to the time scale of the high-velocity 
jet. The third case is the interpretation proposed by Kahane et al. (1996), in which the 
medium- velocity wind came from a bipolar outflow with large opening angle. However, the 
orientation of the approaching and receding parts of the medium- velocity wind suggests that 
the axis of this bipolar fiow should be close to the plane of the disk. If we take into account 
the geometry of the low- velocity disk and the high- velocity jet, the case 3 scenario is less 
hkely, although we cannot rule out the possibihty that V Hya ejects two orthogonal outfiows 
as in the case of CRL 2688 (e.g.. Cox et al. 2000). In Figure 4, we show a schematic picture 
of these three kinematic components based on the scenarios of case 1 or 2. 

The origin of the disk and jet system like V Hya has not yet been elucidated. However, 
theoretical models (e.g., Soker 1997; Soker & Rappaport 2000) suggest that it is difficult for 
a single star to form a flattened disk-like structure. The evidence of rapid rotation of V Hya 
(Barnbaum, Morris, & Kahane 1995) suggests that this star has a companion. It is likely 
that the close companion plays an important role in enhancing the equatorial mass loss, and 
producing the disk and jet structure (e.g., Soker 1992; Morris 1987). The circumstellar 
structure of V Hya resembles that of the proto-PN, CRL 618, which is also associated with 
large equatorial tori and extremely high- velocity bipolar outflow (Sanches Contreras & Sahai 
2004) . The circumstellar structure of V Hya imphes that the asymmetric structure with disk 
and jet appears when the central star is still in the AGB phase. 

We wish to thank all the SMA staff in Hawaii, Cambridge, and Taipei for their enthusi- 
astic help during these observations. We also thank P.T.P. Ho, J.M. Moran, and anonymous 
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Fig. 1. — (a)-(i) Velocity-channel maps of the CO J— 3-2 line. The coordinates of the map 
center are aj2ooo — 10^ 51™ 37.25®, 5j2ooo — -21°15' 00.5". The velocity range of each panel is 
shown in a bracket. Cross denotes the position of the continuum source (see panels k and 1). 
The contours are drawn every 3a. The 3a level of each map is 4.2 Jy beam~^ km s~^ in (a), 
3.3 Jy beam~^ km s~^ in (b) and (h), 2.7 Jy beam~^ km s~^ in (c)-(g), and 4.7 Jy beam^^ 
km s~^ in (i). (j) A map of CO J=3-2 integrated over the velocity range from —44.9 to 
+12.2 km s~^. The contours are drawn every 3a (13.7 Jy beam~^ km s~^) step, (k) A map 
of the 335 GHz continuum emission. The lowest contours are 60 mJy beam~^ {3a) with the 
step of 20 mJy beam~^ (l^")- (1) A map of the 241 GHz continuum emission . The contours 
are every 9 mJy beam~^ (la) with the lowest contour of 27 mJy beam~^ (Scr)- 
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Fig. 2. — Position- velocity map of CO J— 2-1 along the east-west cut through the field center 
in the three different velocity ranges. The effective velocity and spatial resolutions are 10.6 
km s^^ and 5.0", respectively in (a) and (c), and are 2.1 km s""^ and 5.0", respectively in (b). 
The contour interval is 0.055 Jy beam^^ (la) with the lowest contour at 0.11 Jy beam~^ 
(2(t) in (a) and (c). The contours in (b) are drawn in a log scale; 0.40, 0.50, 0.63, 0.80, 1.00, 
1.3, 1.6, 2.0, 2.5, 3.2, 4.0, 10.0, 12.6, 15.9 Jy beam"!. 



- 10 - 



Vlsr 



-40 -30 -20 -10 10 

^' ' ' ' ^ 




20 10 -10 -20 

R.A. Offset (arcsec.) 



Fig. 3. — The CO J=2-l maps integrated over the velocity range of ±(60.0-162.2) km 
from the systemic velocity of —17.5 km s^^ (blue and red contours) superposed on the 
intensity-weighted mean velocity map of the CO J=2-l in the velocity range of ±30 km s~^ 
from the systemic velocity. The contours are drawn every 14.3 Jy beam~^ km s^^ (la) with 
the lowest contours of 42.9 Jy beam~^ km s~^ (3cr). 
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Fig. 4.— 



Schematic picture of the circumstellar structure of V Hya. 



